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ABSTRACT 


The  solution  is  given  for  the  problem  of  screw  dislocation  in 
hexagonal  crystals  with  long  range  interatomic  interactions. 

The  field  equations  of  nonlocal  elastic  solids  are  employed 
to  determine  the  anti-plane  shear  stresses  and  the  elastic 
energy  for  a screw  dislocation  in  the  basal  plane.  Interestingly 
enough,  none  of  the  classical  stress  and  energy  singularities 
are  present  in  the  nonlocal  solutions.  Maximum  shear  stresses 
are  calculated  for  several  hexagonal  crystals  and  compared 
with  the  isotropic  materials.  Theoretical  shear  stress  to 
initiate  a dislocation  with  a Burger's  vector  of  one  atomic 
distance  is  calculated  and  found  to  be  in  the  acceptable  range 
known  from  the  lattice  dynamic  calculations. 


1 . INTRODUCTION 

The  classical  elasticity  solutions  of  Volterra  dislocations  are  well 
documented  (cf.  II)).  The  simplest  among  these  is  the  screw  dislocation 
in  an  Isotropic  elastic  solid  which  exhibits  a 1/r -singularity  for  the  shear 
stress  and  a logarithmic  singularity  for  the  elastic  energy.  The  classical 
elasticity  solutions,  therefore,  fail  in  a "core  region"  near  the  center 
line  of  a cylinder.  In  a previous  paper  [2]  we  have  shown  that  such 
singularities  are  not  present  in  the  solutions  based  on  the  nonlocal  elasticity 
theory.  This,  rather  recent  theory  [3,4],  models  the  solids  much  more 
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satisfactorily  in  that  the  effect  of  long  range  interatomic  interactions  are 
taken  into  account  in  the  stress  constitutive  equations.  This  in  return 
allows  the  treatment  of  geometrical  discontinuities  and  those  associated  with 
the  physical  inputs  (force  distributions,  wave  lengths,  energy,  etc.)  in  a 
more  satisfactory  manner.  Yet,  it  is  a continuum  theory  so  that  all  problems 
can  be  formulated  as  boundary-value  problems.  The  failure  of  classical 
elasticity  theory  in  the  dislocation  core  region  has  led  physicists  to  invent 
various  atomistic  models  to  provide  an  estimate  for  the  state  of  stress  and 
energy  in  this  region. 

Encouraged  with  the  results  of  our  recent  work  on  isotropic  nonlocal 
elastic  solids  [2]  we  felt  that  the  analysis  should  be  extended  to  anisotropic 
solids,  since  in  anisotropic  materials  the  solution  is  influenced  highly  with 
the  orientational  effects.  Hexagonal  crystals  represent  technically  an 
important  class  of  materials  with  well-known  stable  dislocation  patterns. 

Here  we  give  the  solution  of  the  screw  dislocation  problem  with  the  Burger's 
vector  lying  in  the  basal  plane.  In  section  2 we  present  a brief  summary  of 
the  field  equations  of  the  nonlocal  elasticity  theory.  In  section  3 we 
obtain  the  solution  of  the  screw  dislocation  problem  leading  to  expressions 
of  the  stress  fields  and  the  elastic  energy.  In  section  4 some  results  of 
computer  calculations  are  presented  and  the  maximum  shear  stresses  that  cause 
a single  dislocation  in  several  hexagonal  crystals  (Mg,  Apatite,  Cd , 7.n ) 
are  calculated.  The  shear  stress  distributions  along  a radial  plane  for 
these  crystals  display  considerable  differences  from  those  for  the  isotropic 
solids.  Crat ifyingly  no  stress  and  energy  singularity  occur  so  that  by  use 
of  the  maximum  shear  stress  hypothesis  it  is  possible  to  calculate  theoretical 
shear  stress  (cohesive  stress). 
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2 . FORMULATION 

The  basic  equations  of  linear,  homogeneous,  nonlocal  elastic 
solids,  in  the  static  case  with  vanishing  body  force,  are  [3,4]: 


(2.1) 


0 


(2.2) 


'kl 


akimn(x,’x)emn(x')dv(^)  ’ 

1/ 


(2.3) 


mn 


= i(u  +u  ) 
m,n  n,m 


where  the  only  difference  from  the  classical  elasticity  theory  is  in  the 
stress  constitutive  equations  (2.2)  which  state  that  the  stress  at  a 
reference  point  x is  a function  of  strains  at  al 1 points  x'.  As  usual  we 
use  the  summation  convention  for  the  repeated  indices  over  the  range  (1,2,3) 
and  denote  the  partial  differentiation  with  respect  to  x^  by  a comma,  e.g. 

‘kl.m  E 3tk  */9x™  ’ uk,*(x,)  = 3uk/3x,H  ’ 


The  total  strain  energy  of  the  body  is  given  by 


(2.4) 


T.  = i 


tkielkdv 


akemn(x,*x)ek^(x)emn(x,)dv(x’)dv(x) 


1/  V 


From  this  expression  and  (2.2)  it  is  clear  that  the  nonlocal  elastic  moduli 
aktmn^x',}0  possess  the  symmetry  regulations 


(2.5) 


akfmn(x’*x) 


a„.  (x',x) 

Hkmn  _ 


a , „ (x,x' ) 
rankl 


anmk5.(x’X,) 
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Hence  there  are  21  independent  functions  for  the  nonlocal  anisotropic  solids 
that  contribute  to  the  total  elastic  energy.  The  translational  invariance 
of  (2.2)  dictates  that  for  homogeneous  materials  a must  depend  on  x'-x. 


(2.6) 


a,  (x' ,x) 
kcmn 


akemn[(xl 


)2  , (x2~x2)2  , (x'}-x.3)2  ] 


In  conformity  with  the  phonon  dispersion  curves  in  lattice  dynamics  we 

can  approximate  a.  , (x'-x)  by 

kkmn  . _ 7 


(2.7) 


a.  (x'-x) 
k.  mn 


Ckfmna|(xl 


-Xj ) ",  (x2— x2)2  , (x^-X.j)  2 ] 


where  art’  the  elastic  constants  of  the  classical  elasticity  and 

i(x'-x)  is  an  attenuation  funct^ioji  that  depends  on  x'-x. 

From  the  physics  of  solids  the  following  properties  of  a(x'-x)  are 
obv ions. 

(a)  i(x'-x)  attenuates  rapidly  with  x^-x^. 

(b)  In  the  classical  elasticity  limit  a(x'-x)  must  become  a Dirac 
delta  measure. 

Based  on  these  observations  we  assume  that: 

(i)  t(x'-x)  is  a continuous  function  of  x'-x,  with  a bounded  support 

fl  where  a>0. 


(ii) 


(2.8) 


’ 

u(x' -x)dv(x' ) 

U 


1 . 


In  an  exactly  similar  fashion  to  our  work  [2]  we  can  now  prove  that  (2.1) 
and  (2.2)  are  satisfied  if  and  only  if 
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(2.9) 


o.  , = 0 
k tf , k 


in  V , 


where 


(2.10) 


ki! 


c,  „ e 
k f.mn  mn 


is  the  classical  Hookes'  law.  This  result  may  be  stated  as: 

Theorem  of_Co r respondence : 

The  displacement  field  of  the  boundary-value  problem  of  nonlocal,  anisotropic 
elasticity,  under  the  assumptions  (i)  and  (ii),  is  identical  to  that  of  the 
classical  anisotropic  elasticity  theory.  The  stress  field  is,  however, 
given  by 


(2.11) 


Ckf(x) 


j a(x'-x)akf(x') dv(x')  . 

V 


These  results  are  valid  for  any  kind  of  anisotropic  solid. 

For  the  hexagonal  crystals  the  number  of  non-zero  independent  elastic 

constants  c,  is  five  (cf.  [1]  p.  428)  so  that  (2.10)  has  the  form: 
kP.mn  1 ' 


°11 

C11 

C12 

C13 

0 

0 

0 “ 

" en“ 

°22 

r12 

C22 

C12 

0 

0 

0 

e22 

°33 

z 

C13 

C12 

C11 

0 

0 

0 

e33 

CM 

C 

0 

0 

0 

r44 

0 

0 

e23 

°31 

0 

0 

0 

0 

c55 

0 

e31 

°1 2 

0 

0 

0 

0 

0 

n 

i 

e12 

r 
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where 


C11  C11 11  ’ C 12  “ C 1 1 2 2 ’ C13  = C 1 1 33  ’ C22  C2222  , 

(2.13) 

C44  = c2323  ’ c 55  = i(cH_c13) 


The  attenuation  function  a (x'-x)  may  also  have  different  rates  of  decay  in 
the  basal  plane  perpendicular  to  the  x^-axis  than  in  the  x^ -direct  ion . 
Considering  the  fact  that  in  the  basal  plane,  hexagonal  crystals  are 
isotropic  and  recalling  (2.7)  we  must  have 


(2.14)  a = a [ (x'-x  ) (x'-x  ) , (x^-x2)  ] , 6=1,3  . 


The  specific  form  of  this  function  can  be  fixed  by  the  dispersion  curves 
available  in  lattice  dynamics.  A very  useful  one  is  the  function' 


(2.15)  a = aoexp[-(k] /a)2(xj^-x(?)  (Xp-Xg  )-(k2/a)2(x^-x2)^]  ^ B=1  3 


where  a.  is  determined  by  the  normalization  (2.8).  The  constants  k^  and 
k,,  are,  respectively,  the  attenuation  factors  in  the  basal  plane  and  along 
the  x^-axis  and  a is  the  lattice  parameter. 


i-' 


While  this  function  does  not  have  finite  support  nevertheless  it  gives 
the  Dirac  delta  measure  in  the  limit  a-*0. 
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3.  SCREW  DISLOCATION 

A pure  screw  dislocation  in  the  z-direction  in  the  basal  plane  of  a 
hexagonal  crystal  is  possible.  Fig.  1.  Referred  to  rectangular  coordinates 
(Xj=x,  x.,=y,  x^=z),  the  nonzero  components  of  the  displacement  and  stress 
fields  are  given  by  (cf.  [1),  p.  426) 

(1.1)  u^  = ~ tan*  (B  ^v/x) , 


(3.2)  a 


xz 


Ab 

2n 


X 

2 2 ’ 
Bx  +y 


i 1.3) 


o = ABb  x 
vz 


2tt  o 2^  2 
Bx  +y 


where  b=(0,0,b)  is  the  Burger's  vector,  and  A and  B are  constants  related  to 


the  classical  elastic  moduli  c„  by 


(3.4) 


A = 


:c44(circi3)/21' 


B 2c44/(circ13) 


The  elastic  energy  per  unit  length  of  a cylinder,  with  inner  and  outer  radii 
rQ  and  R,  is  given  by 

2 

(3.5)  1/ L = In  (i) 

4tt  r 

o 


From  these  results  it  is  clear  that  the  stress  field  has  a 1/r-singularity 
and  the  energy  has  a logarithmic  singularity  as  r^O.  This,  of  course,  is 
the  troublesome  state,  well-known  in  classical  elasticity. 

According  to  the  theorem  of  correspondence  we  can  employ  these  results 
to  obtain  the  solution  in  nonlocal  elasticity.  To  this  end  we  first  caleulat 
the  classical  stress  fields  in  cylindrical  coordinates.  Thus, 
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(3.6) 


...  Ah  1 

o . = -a  sin0+o  cos0  = x— 

0z  xz  yz  2tt  r 

A (B-l)b  1 tan0 
rz  2n  r B+tanTe 


whore 


cosO  = x/Cx'+v")5  , sin0  = y/(x^+y‘)^ 


In  curvilinear  coordinates  the  physical  components  of  the  stress  field 

(k) 

t ^ j is  given  by  (see  [2]) 


(3.7) 


J u(x'-x)0(k  )(f^(x')6t'  i6kk,dv(x') 


i ' k. 

where  5 and  <5  . * are  the  direction  cosines  between  the  curvilinear 
l k 

k k 

coordinates  x'  and  x . In  this  case  we  have 


« ' 1 . = «J2,  = er-e’r  = e0-e’0  = cos(0'-0)  , = e^e’,  = 1 


(3.8) 


*5"  j i = ee’e'r  = sin(e'-O)  , 6 1 , i = e r • e g = -sin(0'-0). 


where  (er,  e()  , e^)  are  the  unit  vectors  of  the  cylindrical  coordinates  at  x 
and  ( t*  * r , e'0>  e'z)  are  those  at  x'.  Fig.  2. 

Substituting  (3.8)  into  (3.7)  we  obtain 


(3.9) 


t rz  e j a(x'-x)  [ -o0z (x' )sin(0 '-0)+orz (x ' )cos(0 '-0) ]dv(x ' ) , 


t = 'i(x'-x)  [a.  (x ' ) cos (0 ' -0 )+o  (x ' )sin(0 ' -0) ]dv(x ' ) 

Hz  I _ _ 0z  _ rz. 
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ilio  normalization  constant  < in  (2.  15)  is  obtain.  .1  ' . inr.  ■ 
over  the  infinite  space.  This  gives 


( 1. 10) 


3 


In  cylindrical  coordinates  for  the  function  a we  have 

«(x'-x)  = exp[-(k1 /a)"(z'-z)2l •exp[-(k) /a)'  r“cos~n 

2 2 2 o?9 

-(k2/a)  r sin  0]  *exp[-(kj/a)  (r1  cos  f) 1 -2rr 'cos  jcoso  1 ) 

-(k9/a)^(r,,'sin‘'0'-2rr'sin0sin0')] 

It  we  note  that  the  function  a is  even  in  ( 0 ' — Q ) the  integrals  in  0.9) 
take  the  forms 


= A(B-Pb 

rz  2 it 


r°° 

0 J 

tan6 


a(x'-x)  — cos(e'-e)d6'dr'dz’ 

3 ' ' B+tan  9’ 


(3.12) 


Sz  2rr 


Ab 

r 

2^  j 

J 

L 

2it  ( 

*(x'  -x ) | cos < ' -••)  + ( B- 1 ) — / sin (6-  Olde’dr’dz* 

o B+taii  0' 


Integration  in  (3.12)  over  r'  and  z'  can  be  carried  out  leading  to 


Akir 

rz  = 3/2 

sir 


2 TT 

r / r.  sin9'cos0'  , , 

j(B-l)  cosO’-  'U  (,  ,B,e')d6 

o 1+(B-1 )cos  0 ' 


(3. 13) 


"0Z 


4 TT 


AkK 

3/2 


2 TT 


•os  (n '-fi)  + ( B- 1 ) jAne'cose'  >ln(e,  e)]f((  . 

l+(B-l)cos  9' 
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, k=kj  , k ^ /k ! = 

exp(-  “|  l+(*f ^-llsin^Q 


2 2 
(cos  _+>  iin8sln  ’ ) ' 

2 2 

1 + (•"“- 1 ) s in “ 0 ' 


, , cosBcosB '+►'  sinOsinB'  .. 
1+erf  [<  , - ]} 

[1+(k  -l)sin  0 ' ] ■’ 


( l + (tc2-l)sin20'  I * 


• :i  ;i  I .i.  i.*ment  I ield  is  given  hv 


• !.!'>)  u = ti  =0  , u = ■=—  tan  (B  ‘tanO) 

r H z 2tt 


so  that  the  onlv  non-zero  componenet  of  the  strain  tensor  is 


n’lh)  en  = ^ |l  + (B-l)cos20l"'  . 

•z  s 7t  r 


Tin  total  strain  energy  per  unit  length  in  the  z-direction  is 
calculated  bv 


(3.17) 


[ 1+(B-1 )cos‘ 0 ] *t  (r,0)drd0 

0Z 


4.  DISCUSSION 

The  shear  stress  t takes  its  maximum  values  at  The  non- 

0Z 

dimensional  shear  stress. 


(4.1) 


T - t /t  t = ~ - 

0z  C0Z/  o ’ 0 ■ 2tt  a ’ 


1 1 


is 

We 


p 1 ot  t t'd  in 
have  select 


Fig.  1 as  a function  of  p for  various  hexagonal  crystals, 
ed  the  value  of  < as  the  ratio  of  elastic  constants  cr))/c 


1 1 


with  the  consideration  that  the  attenuation  in  a given  direction  is  probably 

proportional  to  the  elastic  modulus  in  that  direction.  The  elastic 

constants  used  are  taken  from  [5]  and  they  are  listed  in  Table  1 together 

with  > and  the  maximum  value  of  t and  the  p=p  at  which  it  occurs.  The 

Oz  m 

t 't^  is  also  shown  on  Fig.  3 as  a function  of  p for  the  isotropic  crystals. 

it  is  clear  that  the  maximum  shear  stress  and  its  location  is  greatly 

affected  by  the  anisotropy.  In  particular,  for  Cd  and  ?.n  the  maximum 

shears  are  nearly  one  half  of  that  for  the  isotropic  solids.  In  Figs.  4 

nd  5 we  display  the  non-dimensional  shear  stresses  t.  and  t =t  / t as 

9z  rz  rz  o 

functions  of  9 for  magnesium.  The  maxima  occurs  at  p-1.10,  9=it/2  for 

t and  p-1.40,  9=m/4  for  t 
z rz 

For  engineering  purposes  it  may  be  useful  to  give  the  ratios  of  shear 

stresses  to  those  for  the  isotropic  solids.  These  are  shown  in  Fig.  6 as 

fun  t i ms  of  p at  6=m/2.  These  curves  may  be  used  to  make  estimates  for 

ether  hexagonal  crystals  with  different  x.  Since  t is  much  smaller  than 
* r z 

t we  have  not  provided  corresponding  curves  for  t^. 

Unlike  the  results  in  the  classical  theory  the  shear  stresses  possess 

no  singularity  but  acquire  maxima.  Consequently  we  can  equate  the  maximum 

shear  stress  to  the  cohesive  shear  stress  to  obtain  the  condition  to  produce 

a dislocation  of  single  atomic  distance.  Thus  taking  b/a=l  we  have 

calculated  the  ratio  of  theoretical  shear  stress  to  the  shear  modulus  c(/. 

44 

These  are  listed  in  Table  2 for  different  attenuation  constants  k=l , 1.25 
and  1.50.  The  value  k= 1 . 50  makes  the  dispersion  curve  for  the  plane  shear 


waves  obtained  theoretically  by  using  (2.15)  nearly  coincident  with  those 
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obtained  from  experiments  of  Joynson[6],  It  appears  that  the  theoretical 
shear  stress  calculated  is  about  twice  the  value  based  on  lattice  dynamics 
calculations  (cf.  [7J^  p.  19)  for  X.n.  Considering  the  inaccuracies 
involved  in  the  estimate  of  interatomic  force  laws,  it  is  clear  that  the 
present  results  are  in  the  right  range.  In  fact,  if  one  takes  k somewhat 
smaller  than  1 and  k slightly  different,  it  is  possible  to  lower  these 
ratios.  However,  one  cannot  place  any  great  faith  in  these  values  in  the 
absence  of  experiments  sufficiently  accurate  for  the  atomic  scale  phenomena. 
The  total  energy  given  by  (3.17)  may  be  expressed  as 

(4.2)  E = Ab2L£ 

o 

where  E depends  on  «,B  and  P=kR/a.  For  various  materials  E is  calculated 
o o 

and  listed  in  Table  3.  The  case  of  isotropic  materials  agree  very  well 
with  the  result  given  in  [2]^.  Er  grows  indefinitely  with  the  radius  R 
becoming  infinite  for  R=oo  as  expected. 


'in  [2]  eqs.  (4.2)  and  (4.3)  contain  a typographical  error.  Right-hand 
side  of  E should  be  multiplied  by  i in  both  (4.2)  and  (4.3). 


Table  1.  Maximum  Shear  Stress 


1 3 


Elastic  Constants 

Material 

C11  } c 1 2 | c13  | 1 22  | c44 

K = 

T 

c - 

m 

xlO*  *dyn/cm^ 

(c22/cll)  J 

(t  /t  ) 
6 z o 

(kr  /a) 

L--v  - 

Zn 

16.5 

5.0 

3,1 

6.2 

- 

3.96 

0.376 

| 

0.3370 

2.65 

Mg 

5.93 

2.14 

2.57 

6. 15 

1 . 64 

1.037 

0.6555 

1.08 

Cd 

11.4 

4.0 

3.94 

5.08 

1 

2.0 

0.446 

0.3911 

2.25 

Apat i te 

16.7 

6.6 

1.31 

14.0 

6.63 

0.838 

0.5871 

1.29 

Ice 

1.34 

J 

0.53 

1 

0.65 

1.45 

0.313 

1.082 

L 

0.6824  | 

1.03 

I 


I 

I 


Table  2.  t /c.. 

0z  44 


k 

1.00 

1.25 

1.50 

Zn 

0.0698 

0.0872 

0. 1049 

Mg 

0. 1056 

0. 1320 

0. 1584 

Cd 

0.0850 

0. 1062 

0.1275 

Apat ite 

0.1007 

0. 1258 

0.1510 

Ice 

0.1140 



0.1425 

0. 1710 

Table  3:  Coefficient  E of  Total  F.nergv;  E =E/Ab 
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